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ABSTRACT 



Aims. We compute and analyze the evolution of primordial stars of masses at the ZAMS between 5 Mq and 10 Mq , with and 
without overshooting. Our main goals are to determine the nature of the remnants of massive intermediate-mass primordial 
stars and to check the influence of overshooting in their evolution. 

Methods. Our calculations cover stellar evolution from the main sequence phase until the formation of the degenerate cores and 
' the thermally pulsing phase. 

f***) , Results. We have obtained the values for the limiting masses of Population III progenitor stars leading to carbon-oxygen and 
oxygen-neon compact cores. Moreover, we have also obtained the limiting mass for which isolated primordial stars would lead to 
core— collapse supernovae after the end of the main central burning phases. Considering a moderate amount of overshooting the 
mass thresholds at the ZAMS for the formation of carbon-oxygen and oxygen-neon degenerate cores shifts to smaller values by 
Q ' about 2Mq. As a by-product of our calculations, we have also obtained the structure and composition profiles of the resulting 
S—j \ compact remnants. 

CZ) . Conclusions. Opposite to what happens with solar metallicity objects, the final fate of primordial stars is not straightforward 
' determined from the mass of the compact cores at the end of carbon burning. Instead, the small mass-loss rates typically asso- 
ciated to stellar winds of low metallicity stars might allow the growth of the resulting degenerate cores up to the Chandrasekhar 
V~J ■ mass, on time scales one or two orders of magnitude shorter than the time required to loose the envelope. This would lead to 
' the formation of supernovae for initial masses as small as ~ 5Mq. 
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1. Introduction 

Primordial stars — that is, stars formed in the early 
Universe — are direct inheritors of the matter synthe- 
sized during the Big Bang (Alpher & Herman, 1950; Olive, 
2000). The initial composition of these objects — which 
is characterized by the absence of metals — determines 
their evolution and, ultimately, the fundamental proper- 
ties of the resulting stellar remnants, as well as the amount 
and composition of the matter returned to the interstellar 
medium after the main evolutionary phases. Therefore, 
the evolution of these stars is important for a correct 
modeling of the population of the degenerate remnants of 
Population III stars and, hence, for a better understanding 
of the distribution of baryonic dark matter, whereas the 
significance of the ejected material rests on its influence 
on Galactic chemical evolution. 
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Even though the detection of primordial composition 
objects — that is stars with [Fe/H] < — 8.3 — has not 
been possible up to now, during the last few years several 
encouraging observations have been made. For instance, 
Bessel & Norris (1984) detected a red giant with [Fe/H] ps 
—4.6 and Christlieb et al. (2002) measured the metallicity 
of HE 0107-5240, a low mass star with [Fc/H] w -5.3. 
Even more recently, Frebel et al. (2005) have observed 
HE 1327-2326, an unevolved low-mass object, which was 
found to have [Fe/H] w —5.4. Consequently, both the in- 
trinsic theoretical interest and the growing observational 
evidence for the existence of primordial stars make the 
study of these objects one of the hottest topics in stellar 
astrophysics nowadays. 

The exact shape of the initial mass function of 
Population III stars is far from being known. The stan- 
dard view advocates for a strong bias toward very massive 
stars, say between 100 and 600 Mq — see, for example, 
Bromm, Coppi & Larson (2002). Nevertheless, arguments 
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have been put forth recently that point to the possible ex- 
istence of a population of very metal poor stars of low- and 
intermediate-mass. For instance, from the theoretical side, 
Nakamura & Umemura (1999, 2001) have simulated the 
growth of instabilities in protostellar primordial clouds in 
2 and 3 dimensions. Their results indicate that the initial 
mass function could be bi-modal with a peak at ~ 100 Mq 
and a secondary peak at a few solar masses. Additionally, 
Machida et al. (2005) have shown that the explosion of 
primordial supernovae can drive the Jeans masses of pri- 
mordial clouds from a few 10 6 Mq to values well below 
1 Mq, hence allowing for the formation of intermediate- 
mass stars, but with non-negligible metallicities (of the 
order of [Fe/H] sa —3.0). Furthermore, Johnson & Bromm 
(2006) have recently raised more theoretical considera- 
tions that suggest that the formation of primordial low- 
and intermediate-mass stars is viable. On the other hand, 
and from the observational side, Schneider et al. (2006), 
on the basis of three tentative NICMOS UDFs sources at 
z s=s 10 and the derived WMAP electron scattering optical 
depth, predict a rather standard form of the IMF, with 
stars within the 1 Mq ^ M ^ 100 Mq range. 

The evolution of zero metallicity stars has been ana- 
lyzed in several recent publications. Just to mention some 
of them, and without trying to be complete, it is important 
to realize that Heger et al. (2001) and Heger & Woosley 
(2000) have focused on massive objects, whereas Limongi 
et al. (2000), Marigo et al. (2001), Chieffi et al. (2001) 
and Siess et al. (2002) have thoroughly studied the evo- 
lution of low- and intermediate-mass stars. However, the 
evolution of massive intermediate-mass primordial stars 
has been largely ignored, except for the case of the evo- 
lutionary sequence of a 9 Mq star presented in Gil-Pons 
et al. (2005). Probably, one of the reasons for this is the 
heavy computational load involved in computing the evo- 
lution during the carbon burning phase for this mass range 
(Ritossa et al. 1995; Garcfa-Berro et al. 1997; Iben et al. 
1997; Ritossa et al. 1999). This is in sharp contrast with 
the situation for stars of Population I and II, that have 
been thoroughly analyzed. For instance, the mass thresh- 
olds for objects that lead to different types of remnants 
have been determined for Population I and II stars by 
Heger et al. (2003), as well as by Eldridge & Tout (2004), 
who computed a rich grid of stellar models ranging both 
in mass — from 5 to 200 Mq — and metallicity — from 
Z = 10~ 5 up to Z = 0.05. The latter authors also payed 
special attention to the adopted mass-loss rates, including 
their variation with the metallicity, and to the effects of 
overshooting. Finally, it is worth mentioning as well that 
the influence of Population III stars on the formation and 
ejection of isotopes to the intergalactic medium has been 
extensively studied by Goriely & Siess (2002) and by Abia 
et al. (2001). 

In this paper, we describe the final stages of the 
evolution of massive intermediate mass primordial stars. 
Reliable results on this matter are essential inputs for 
Galactic chemical evolution models, as well as for the 
supernova theory. Particularly, thermonuclear supernovae 



(SNela), in which a carbon-oxygen white dwarf overcomes 
the Chandrasckhar limit, must rely on accurate initial— to- 
final mass relations to reproduce the observed supernova 
rate, and its dependence on metallicity could represent 
a source of diversity on classical SNela at different red- 
shifts. Specifically, we outline the evolution of the pro- 
genitor stars and determine the mass thresholds to ob- 
tain carbon-oxygen (CO), oxygen-neon (ONe) degener- 
ate cores and core collapse supernovae at the end of the 
main central burning stages. The stars hosting CO degen- 
erate cores could give raise to CO novae or to SNela if 
belonging to a close binary system; those hosting ONe de- 
generate cores could also form ONc novae or supernovae 
by the accretion-induced collapse mechanism (Canal & 
Schatzman 1976; Gutierrez et al. 1996; Gutierrez, Canal & 
Garci'a-Berro, 2005), and core-collapse supernovae would 
directly form from isolated stars when nuclear burning 
proceeds all the way to the formation of an iron peak ele- 
ment dominated core. 

An important part of our results consists in justifying 
that the formation of supernovae docs not require ZAMS 
masses as high as in the case of solar metallicity objects 
but, instead, 6 Mq stars (and probably even less mas- 
sive ones) are enough to produce a supernova explosion. 
This possibility had already been suggested, for instance, 
by Arnett (1969), by Iben & Renzini (1983), by Willson 
(2000) and, more recently, by Zijlstra (2004). The reason 
for such peculiar behavior has to do with the character- 
istics of stellar winds in primordial stars. Low metallicity 
stars are supposed to drive weaker stellar winds (Bowen 
& Willson 1991; Willson 2000) than their solar metal- 
licity counterparts. However, some caution is mandatory 
here, as the cited works are theoretical with phenomeno- 
logical character, and observations do not fully support 

nor discard — their findings yet. Therefore, the de- 
generate cores are able to grow significantly due to shell 
burning before the envelope is lost. As we will show be- 
low, the cores resulting from the evolution of massive 
intermediate-mass primordial stars will be able to grow 
up to the Chandrasekhar mass and, therefore, undergo a 
supernova explosion. 

Of course, all these results depend sensitively on the 
adopted prescription for convection and, in particular, on 
the degree of overshooting. Overshooting represents one of 
the main unknowns in stellar evolution, as observations in- 
dicate more mixing than predicted by standard, local and 
non-local convective prescriptions. Rotation and purely 
three-dimensional hydrodynamic phenomena could ex- 
plain this extra mixing, but this is still the subject of an 
active, ongoing debate. Hence, we have performed calcula- 
tions both with and without overshooting in order to clar- 
ify its effects on the initial-to-final mass relation. Previous 
works on overshooting — see, for example, Young et al. 
(2003) and Demarque et al. (2004) for recent, enlightening 
studies — indicate that overshooting produces an increase 
in the mass of the nuclear exhausted core, a change of its 
internal composition, and a modified chemical composi- 
tion of the surface layers. All these effects are relevant 
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Fig. 1. Evolution in the logp c — logT c plane of our model 
stars up to the point where carbon is ignited. The solid 
lines correspond to the case in which no overshooting was 
taken into account, whereas the dotted lines correspond to 
the calculations in which overshooting was incorporated. 
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Fig. 2. Evolutionary tracks in the Hertzsprung-Russell di- 
agram of our model stars. The solid lines correspond to 
the case in which overshooting was disregarded, whereas 
the dotted lines show the evolution for the case in which 
overshooting was due taken into account. 



for the purposes of our work. In particular it is of the 
maximum importance the extra mixing in the surface lay- 
ers, as population III stars do not experience a significant 
mass-loss due to the small opacity of their outer layers. 
The extra mixing due to overshooting could increase the 
mctallicity of the outer layers and, consequently, increase 
the mass-loss rates, thus changing the final mass of the 
remnants. 

The paper is organized as follows. In Sect. 2 we 
describe the overall evolution prior to carbon burning. 
Section 3 is devoted to study in depth the different car- 
bon burning episodes. In §4 we comment on the thermally 
pulsing phase for this mass range, whereas in §5 we de- 
scribe the properties of the cores of these stars. Finally, in 
section 6 we elaborate on the possible evolutionary out- 
comes that can be envisaged for these stars, discussing in 
depth the effects of the adopted mass-loss rates, whilst in 
section 7 we summarize our major findings, we draw our 
conclusions and we discuss their significance. 

2. Overall evolution prior to carbon burning 

We have followed the evolution of zero-metallicity mod- 
els with masses at the ZAMS between 5 and IOMq, 
and initial values for hydrogen and helium abundances of 
X = 0.77 and Y = 0.23 respectively. Our calculations have 
been performed both taking into account and neglecting 
overshooting. In doing so, we have used the evolution- 
ary code previously described in Gil-Pons et al. (2005), 
which is, itself, an update of the code described in detail 
in Ritossa, Garcfa-Berro & Iben (1995). We remark that 
the opacities are obtained by interpolating in the tables 



of Iglesias & Rogers (1993). The only novelty with re- 
spect to the previous implementations of our code is that 
for the calculations reported here we have incorporated 
overshooting. This will allow us to isolate the differences 
introduced by the treatment of convective boundaries in 
the evolutionary calculations presented here. Specifically, 
the prescription used for the calculation of the convective 
edges is the one described in Eldridge & Tout (2004), in 
which the standard Schwarzschild criterion is modified to 
let convection be the main energy transport mechanism 
in those zones of the star where V ra d > V a d — 5, with 



2.5 + 20C+16C 2 

where £ = P ra d/-P g as and <5ov is set to be equal to 0.12 
(Schroder, Pols & Eggleton 1997). 

The evolution prior to carbon burning follows very 
closely the main trends previously described in Gil-Pons 
et al. (2005). That is, core hydrogen burning (CHB) be- 
gins through the pp-chains until small amounts of carbon 
(X(C) = 10~ 10 by mass) form and allow the onset of the 
CNO cycle. The rest of the CHB phase proceeds therefore 
at the higher central temperatures that characterize this 
cycle. Once CHB has ended, core contraction proceeds and 
the conditions for core helium burning are reached in the 
central regions of the star. Those for hydrogen burning 
are attained at the layers just above the hydrogen ex- 
hausted core. Opposite to what happens with solar metal- 
licity stars, all these phases take place at the blue region of 
the Hertzsprung-Russell diagram. Furthermore, hydrogen 
shell burning in a metal-free environment is not able to 
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Table 1. Main nuclear burning times and core masses. 
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Fig. 3. Upper panel: density profile of the degenerate 
carbon-oxygen cores for the 8 Mq model star just before 
C ignition for the case computed without overshooting 
(solid line) and with overshooting (dotted line). Lower 
panel: temperature profiles of the carbon-oxygen cores 
computed without overshooting (solid line) and with over- 
shooting (dotted line) for the same model star. 



produce the overall expansion of the stellar envelope that 
characterizes the red giant branch. 

As mentioned in Gil-Pons et al. (2005), the main differ- 
ences between the evolutionary tracks for Z = and those 
of Z = Zq are the longer time spent during the main 
sequence phase, the higher surface luminosities reached 
during CHB and the larger helium-exhausted core of the 
Z = models. Table 1 shows the duration oh the CHB 
phase (£chb)) the mass of the resulting helium core (Mn e ), 
the time at which carbon burning begins (£cb), and the 
mass of the corresponding carbon-oxygen core. As can be 
seen, the calculations with overshooting yield time scales 
for core hydrogen burning that are about a 12% longer 
when compared to the cases computed without taking into 
account overshooting. The helium cores after the CHB 
phase has been completed are also about 50% larger and 
this is the reason why the core helium burning (CHeB) 
phase lasts about 6% longer in the cases in which over- 
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Fig. 4. Composition profiles of the degenerate CO cores 
of the 8 Mq model star computed with overshooting (thin 
lines) and disregarding overshooting (thick lines) just be- 
fore carbon burning sets in. 

shooting was incorporated. Finally the masses of the CO 
cores (Mco) prior to carbon burning are considerably 
larger as well, that is, about 25 % more massive for the 
cases computed with overshooting. 

Figure 1 shows the evolution of the central tempera- 
ture versus the central density for both sets of models. The 
solid lines correspond to the cases in which overshooting 
was disregarded and the dotted lines to the cases in which 
overshooting was taken into account. For the sake of con- 
ciseness we only show the evolutionary sequences of the 7, 
8 and 10 Mq model stars. As can be seen, only small differ- 
ences appear between both sets of models during the CHB 
phase. However, the core contraction phases that follow 
the CHB and the CHeB phases take place at higher tem- 
peratures and reach lower values for the maximum central 
density in the cases computed with overshooting. 

The evolution in the Hcrtzsprung-Russell diagram of 
the same model stars is shown in figure 2. In the upper 
panel the track of the 7Mq star is shown, and in the 
middle and lower panel we represent the evolution of the 
8 Mq and 10 Mq models, respectively. The cases with no 
overshooting are represented as solid lines and the cases 
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Fig. 5. Temporal evolution of the main structural parameters during the bulk of the carbon burning phase for the 
8Mq model computed without overshooting (upper panel) and with it (bottom panel). The shaded areas represent 
the evolution of the convective regions, the thick solid line depicts the evolution of the carbon burning luminosity and 
the thin solid line shows the helium burning luminosity. 



with overshooting are represented using dotted lines. As 
can be seen, CHB takes place at higher effective temper- 
atures for the models computed with overshooting, and 
the Hertzsprung gap occurs at larger luminosities. In sum- 
mary, the models computed with overshooting behave as if 
they were more massive than their counterparts computed 
without overshooting. 

This can be better seen by comparing the masses of 
the CO cores just before the carbon burning phase (Tabic 
1). Also, the density and temperature profiles just before 
carbon ignition — which are shown in Fig. 3 for the 8 Mq 
model star — help in displaying this distinctive behav- 
ior. Particularly, the model computed with overshooting 
presents a smaller central density and a shallower den- 
sity profile — see the top panel of Fig. 3 — leading to a 
larger region of high densities and, therefore, to a larger 
core than the model computed without taking into ac- 
count overshooting. Note, however, that the central tem- 
peratures (see the bottom panel of Fig. 3) of both models 
are remarkably similar, although the mass of the core is 
rather different. Figure 4 shows the chemical profiles of the 
same model star. The thick lines correspond to the model 
computed with no overshooting and the thin lines corre- 
spond to the model computed when due account is taken 
of overshooting. The former appears to be 0.3 Mq larger 
than the latter, but the mass abundances of carbon and 



oxygen in the central regions of the star are practically the 
same in both cases. All these quantities are particularly 
important because, as it will be shown in section 4 below, 
the development of extensive carbon burning and, ulti- 
mately, the formation of an oxygen-neon white dwarf cru- 
cially depend on them. The mass threshold at the ZAMS 
for the formation of supernovae from single Z = stars 
also depends significantly on the adopted prescription for 
determining the edge of the convective cores. 

3. The carbon burning phase 

Figure 5 represents the temporal evolution of the base of 
the convective envelope (BCE) as well as the luminosi- 
ties associated to carbon burning (thick solid line) and 
helium burning (thin solid line) for the 8 Mq model stars 
computed without overshooting (top panel) and with over- 
shooting (bottom panel), which is a representative case. 
The convective regions associated to carbon burning are 
also shown for both models. We remark at this point that 
neither the 6Mq nor the 7 Mq model stars computed 
without overshooting develop carbon burning extensively 
although for the latter case there is partial carbon burn- 
ing, reaching the carbon luminosity a maximum value of 
only L/Lq ~ 10 — and, therefore, the remnant cores are 
composed mainly of carbon and oxygen. In fact, our cal- 
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culations show that the mass threshold at the ZAMS for 
effective carbon burning in the partially degenerate core 
to obtain an ONe remnant is w 7.8 Mq, when no over- 
shooting is considered. This mass threshold is w 6.0 Mq 
when overshooting is taken into account. 

The first important feature to be noticed in Fig. 5 is 
the fast advance of the base of the convective envelope 
that, for both models, occurs shortly after the end of the 
core helium burning phase and the onset of helium burn- 
ing in a shell — at t ~ 9.757 x 10 14 s for the model com- 
puted without overshooting and at t ~ 10.2907 x 10 14 s 
for the model computed with overshooting. This first ad- 
vance inward, the so-called "second" dredge-up episode 
- although, in fact, is the first dredge-up episode (Gil- 
Pons et al. 2005) — is followed by a very slow penetration 
of the BCE that finally allows for a moderate enrichment 
in metals of the hydrogen-rich envelope. Although most 
of the material dredged-up consists of helium which does 
not contribute to the metal content of the convective enve- 
lope some small amounts of metals are indeed dredged-up. 
The metallicity of the envelope depends sensitively on the 
mass of the star. For instance, the metallicity of the enve- 
lope of the 5 Mq model star computed without overshoot- 
ing is Z cnv « 10~ 9 , whereas for the 9 Mq star we obtain 
Z e nv ~ 10 -4 . The models computed with overshooting 
show the same behavior: Z cnv « 10~ 7 for the 5 Mq star 
and Z CDV w 10~ 3 for the 8 Mq. As it will be shown below, 
this enrichment plays a key role in the final fate of mas- 
sive intermediate-mass primordial stars. The reason for 
this behavior can be understood by studying the evolution 
of the 8 Mq model star computed without overshooting, 
which shows a particular feature associated with the inner 
advance of the convective envelope. For this model the first 
carbon flash at t ~ 9.766 x 10 14 s occurs near the center 
of the star, but the associated convective zone extends to 
layers relatively close to the BCE (at M r ~ 0.9 Mq) and, 
thus, allows for a decrease of the degeneracy parameter of 
the outermost layers of the core. Consequently, the den- 
sity and temperature barriers decrease and this results in 
an additional advance of the BCE. This, in turn, causes 
a enrichment in metals of the convective envelope. For 
the model computed with overshooting this feature is ab- 
sent, since carbon ignition occurs at the center, too far 
from the edge of the carbon-oxygen partially degenerate 
core (which is located at Mbce — 2.12 Mq) to have any 
relevant consequence in the behavior of the convective en- 
velope. 

The 8 Mq model computed without overshooting ac- 
tually undergoes 5 significant carbon flashes, the first of 
which is the most powerful one (~ 10 8 Lq). Once the 
overall degeneracy of the core decreases after this flash, 
the following ones proceed with more moderate luminosi- 
ties, as it usually occurs for solar metallicity models of 
equivalent masses during this phase. As it also occurs 
in solar metallicity stars of the same mass, each ther- 
monuclear flash generates an associated convective shell 
and an abrupt decrease of the He luminosity (see Fig. 5). 
This is a consequence of the fast injection of energy in 



the uppermost layers of the partially degenerate core, 
which consequently expand and cool and, therefore, ef- 
fectively switch-off the He burning shell. For the model 
computed taking overshooting into account, carbon igni- 
tion occurs at the center. This is a consequence of the 
smaller central densities and lower temperatures achieved 
before carbon ignition (see Fig. 3), which lead respectively 
to a smaller degeneracy parameter and to smaller neu- 
trino leaks. Consequently, the temperature profile shows 
a central maximum, in contrast with the situation for the 
8Mq model computed without overshooting, for which 
the maximum temperature occurs relatively far from the 
center. Additionally, we do not obtain as many flashes as 
for the model in which overshooting was not taken into ac- 
count. In fact, the first carbon flash has a smaller strength 
{Lq ~ 10 7 Lq). This also stems from the fact that the 
degeneracy in the core is substantially smaller. After the 
sudden injection of nuclear energy, the core of this model 
star expands and the overall degeneracy decreases, leading 
to a more gentle carbon burning phase, in which carbon 
is burnt in almost stationary conditions. Moreover, since 
the outer edge of the central convective carbon burning 
region is located deep in the core, the effect on the He 
burning shell is substantially smaller and, consequently, 
the He burning luminosity decreases less than in the pre- 
vious case. 

4. The thermally pulsing phase 

For the models which fail to ignite carbon under par- 
tially degenerate conditions a thermally pulsing phase en- 
sues soon after helium burning ceases in the central re- 
gions. For these models the helium-burning shell is very 
close to the base of the hydrogen-rich convective enve- 
lope. For instance, for the 7Mq model sequence com- 
puted without overshooting the helium-burning shell is 
only ~ 4.1 x 10~ 4 Mq away from the H-He disconti- 
nuity. Consequently, the temperature at this discontinu- 
ity increases to the extent that hydrogen-burning is re- 
activated and the thermally pulsing AGB phase (TP- 
AGB) phase ensues shortly afterward. For the models 
which do ignite carbon the thermally pulsing Super-AGB 
(TP-SAGB) phase ensues after the carbon burning phase 
has been completed. This, for instance, is the case of 
the 8Mq model star computed without overshooting. 
Again, in this case the helium-burning shell (located at 
M r ~ 1.049875 Mq) is very close to the H-He discontinu- 
ity (located at a mass coordinate M r ~ 1.08169596 Mq) 
and the hydrogen burning shell resurrects as well. Note, 
however, that Yoon et al. (2004) have recently shown that 
other parameters, like rotation, define the stability of shell 
burning. 

In Fig. [6] we display the helium luminosity for the first 
few thermal pulses of some selected models. Note that 
although the time scales are rather different, the ampli- 
tudes of the pulses are similar. The fact that these stars 
develop the TP-(S)AGB phase in a way similar to that 
of solar metallicity models of analogous masses is related 
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Fig. 6. Temporal evolution of -Ljje during the early TP-(S)AGB phase of the 7Mq model computed without over- 
shooting (upper left panel), the 8Mq without overshooting (upper right panel), the 5 Mq model computed with 
overshooting (lower left panel), and the 7 Mq with overshooting (lower right panel). 



to the dredge-up processes that follow the development 
of shell helium burning. The associated metal enrichment 
of the base of the hydrogen-rich envelope, as well as its 
heating due to the proximity to the helium burning shell 
(HcBS), allow, as we have just explained, for a new onset 
of hydrogen-burning through the CNO cycle, that occurs 
as soon as the distance between the HeBS and the base of 
the hydrogen-rich envelope is about 10~ 4 Mq. From this 
time on, hydrogen and helium burning alternate them- 
selves as main energy suppliers. In the next two subsec- 
tions we give some details about the thermal pulses that 
ensue in stars that have developed either a carbon-oxygen 
or an oxygen-neon core. 



4.1. Models with carbon-oxygen cores 

To keep consistency with the existing literature, we will 
refer to these models as TP-AGB models. The upper and 
lower left panels of Fig. [6] show, respectively, the first 
few thermal pulses of the 7 Mq model computed with- 
out overshooting and those corresponding to the 5Mq 



model star computed with overshooting. These are rep- 
resentative examples of massive intermediate-mass pri- 
mordial stars with carbon-oxygen cores. We recall that 
models with masses at the ZAMS smaller than 7.8 Mq 
never ignite carbon if no overshooting is adopted, whereas 
if overshooting is taken into account the mass limit is 
about 6Mq. Since both models fail to ignite carbon, the 
TP-AGB phase starts at the end of the CHeB phase. 
However, it is worth mentioning that before reaching the 
TP-AGB phase, both models have experienced a brief 
(and weak) carbon-burning phase, with associated lumi- 
nosities which never exceed ~ 10 Lq. This phase is, in any 
case, unable to modify substantially the chemical profile 
of the corresponding cores. As it happens for solar metal- 
licity TP-AGB models, each helium flash is accompanied 
by the expansion and subsequent cooling of the neigh- 
boring upper layers and, therefore, by the switch-off of 
the hydrogen burning shell. The interpulse period for the 
7 Mq model computed without overshooting amounts to 
3200 yr, whereas it is reduced to 2700 yr for the 5 Mq 
model computed with overshooting. As it can be seen the 
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Fig. 7. Temporal evolution of the edges of the convective shell and of the base of the convective envelope during the 
8 th and the 9 th pulses of the 7 Mq model star computed without overshooting (upper left panel), the 14 th and the 
15 th pulses of the 8Mq star computed without overshooting (upper right panel), the 18 th and the 19 th pulses of the 
5 Mq model computed with overshooting (lower left panel) , and the 23 th and 24 th pulses of the 7 Mq computed with 
overshooting (lower right panel). In the figure are given the values of AMcs, the maximum mass of the convective shell, 
AA/h, the mass through which the hydrogen profile moves between pulses, AA/dUj the amount of mass dredged-up 
during pulse power-down and AA/eb, the mass of the outer edge of the convective shell and the base of the convective 
envelope. The mass of the hydrogen-exhausted core is shown as a dotted line. 



average values of the helium luminosity are very similar 
(Lfjc ~ 10 4 £q) for both models, but the maximum he- 
lium luminosity is somehow larger for the 7 Mq model 
computed without overshooting. 

The upper left panel of figure [7] represents the evolu- 
tion of the BCE and the pulse-driven convective zones 
associated to the 8 th and the 9 th helium flashes of the 
7 Mq computed without overshooting, occurring at times 
t ~ 12.3446 x 10 14 s and t ~ 12.3455 x 10 14 s, respectively. 
It can be seen that the fast injection of energy associated 
to each helium flash causes the development of an inner 
convective zone, as it occurs in the case of solar metallic- 
ity model stars. Simultaneously, the expansion and cooling 
that occurs in the neighboring layers produce the extinc- 
tion of the hydrogen burning shell and, as a consequence, 
a fast advance inwards of the BCE. By the end of the 
helium flash, the H-burning shell moves outward, increas- 
ing the mass of the He buffer and a new He flash occurs 



once enough matter is accumulated. However, the pene- 
tration of the BCE into the regions that have been pro- 
cessed by previous pulse-driven convective zones is hardly 
noticeable. Hence, mixing of the material of the hydrogen- 
exhausted core with that of the convective envelope does 
not allow for a significant variation of the metallicity of 
the envelope and does not affect the subsequent evolu- 
tion. In fact, the increase in metallicity that this model 
star experiences during each helium flash is only AZ cnv m 
10 -8 . For the sake of completeness, we have included in 
Fig. [7] the values of the maximum mass of the convective 
shell during the 9 th pulse, (AMcs = 1.61 x 1CT 3 M Q ), 
the mass through which the BCE moves between pulses 
(AMh — 6.75 x 10~ 4 Mq), the mass dredged-up after 
the 9 th helium flash (AM DU = 6.17 x 1(T 4 M ), and the 
thickness of the region between the outermost edge of the 
convective shell and the base of the convective envelope 
(AM EB = 1-12 x 1CT 4 M ). 
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Fig. 8. Profiles of the adiabatic gradient (dotted line) and the radiative gradient (solid line) around the 20 th thermal 
pulse of the 5 Mq model computed with overshooting. Panel a shows the situation near the maximum luminosity 
of the pulse (at t = 2.4341617 x 10 15 s). Panel b depicts the same quantities at at t = 2.4341619 x 10 15 s, which 
corresponds to the maximum extent of the inner convective zone. Panel c shows the gradients near the maximum 
advance inwards of the BCE (at t = 2.4341645 x 10 15 s). Finally, panel d shows the gradients at the end of the pulse 
(at t = 2.4341699 x 10 15 s). 



To this regard it is important to mention here that 
the efficiency of the third dredge-up strongly depends, 
not only on the treatment of convection, but also on the 
numerical treatment. Therefore, it is not surprising that 
the different authors that have studied this evolutionary 
phase have obtained different results. For instance, Marigo 
et al. (2001) did not find any evidence of third dredge-up 
in their calculations, although, as these authors point out, 
this effect could be due to the small number of thermal 
pulses followed. Chicffi et al. (2001) treated the convec- 
tive boundaries according to the prescription of Herwig et 
al. (1997) and, moreover, extended the mixing below the 
convective envelope in such a way that the composition 
discontinuity between the helium-rich and the hydrogen 
rich layers was considerably smoothed. Accordingly, these 
authors found a significant third dredge-up after about 
13 thermal pulses. This result is similar to the one ob- 
tained by Siess et al. (2002). Unfortunately, the treatment 
of convection should be validated by calibrating theoret- 
ical models with observations, which is not possible for 



primordial stars. However, it is worth pointing out that 
even authors whose results for solar metallicity stars sup- 
port the existence of an efficient third dredge-up (Karakas 
& Lattanzio 2002), have found some evidences that might 
point to a less efficient third dredgc-up for the case of 
metal-poor stars (Karakas et al. 2006). Finally, it is also 
important to realize that in our calculations we have not 
included the hot-bottom burning. In fact, it can be ex- 
pected that hot-bottom burning would increase the rel- 
ative abundance in the envelope of nitrogen and oxygen 
with respect to carbon, but not that it changes critically 
the total metallicity of the envelope. So, even though we 
have followed a relatively small number of thermal pulses 
(up to 20 at most), we do not expect any drastic effect on 
our results due to extra pollution by hot bottom-burning. 

The lower left panel of figure 7 represents the behav- 
ior of the convective regions of the 5Mq model com- 
puted with overshooting in the period of time elapsed from 
shortly before the beginning of the 15 th pulse until shortly 
after the 16 th pulse is over. In this case, the intcrpulsc in- 
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Fig. 9. Chemical profiles of the main chemical species for the 5 Mq model star during the 20 th thermal pulse. The 
evolutionary times are those of Fig. [8l The profiles are given at the maximum luminosity of the pulse (panel a), at 
the maximum extent of the inner convective zone (panel b), near the maximum advance inwards of the BCE (panel 
c) and at the end of the pulse (panel d) . 



terval amounts to 2700 yr, somewhat shorter than that 
of the 7 Mq model star previously discussed. Moreover, 
the evolution of the base of the hydrogen-rich convective 
envelope and of the convective zones associated to the he- 
lium flashes is totally equivalent to that already done and 
the surface composition of this model docs not change ap- 
preciably. The values of AM CS , AM H , AM DU and AM EB 
are also displayed in the lower left panel of figure [7] 

In order to get a better insight of the evolution dur- 
ing the occurrence of a helium flash, in figure [5] we have 
represented as a dotted line the adiabatic gradient (V a d) 
and as a solid line the radiative gradient (V ra d) in the 
region near the BCE, for some selected times during the 
20 th thermal pulse of the 5.0 Mq model sequence com- 
puted with overshooting. Panel (a) corresponds to time 
t = 2.4341617 x 10 15 s. At this time the helium luminosity 
is Ln e ~ 2.3 x 10 4 Lq and the convective shell is just form- 
ing. Note as well that for this model the hydrogen luminos- 
ity is rather small Ljj = 4.4 Lq and, therefore, the BCE is 
able to begin its advance inwards. In the zones in which the 
radiative gradient is larger than the adiabatic one, accord- 
ing to the (modified) Schwarzschild criterion used for our 



calculations, there is complete mixing of the chemical ele- 
ments. Panel (b) corresponds to time t = 2.4341619 x 10 15 
s, for which the pulse-driven convective region reaches 
its maximum extension (AM = 1.36 x 10 _3 Mq). Panel 
(c) shows Vad and V ra d at time t = 2.4341645 x 10 15 s, 
when the BCE reaches its maximum extent inwards, at 
M r = 0.9397 Mq. Finally, panel (d) corresponds to time 
t = 2.4341699 x 10 15 s and shows the recession of the BCE 
once the HBS is fully operative again (Ljj = 1-2 x 10 4 Lq). 
As can be seen, the total amount of mass dredged-up is 
very small even if overshooting is taken into account. 

Figure [9] shows the chemical profiles of the main nu- 
clear species at the times already selected for panels (a) 
to (d) of figure [H In panel (a) the 4 He buffer formed due 
to hydrogen shell burning during the previous interpulse 
period can be well seen between the mass coordinates 
M r = 0.93900 Mq and M r = 0.93963 Mq. Later, at time 
t = 2.4341619 x 10 15 s, corresponding to panel (b), the ef- 
fects in the chemical profile of the convective HeBS can be 
seen. Note the very small thickness of the region in which 
He has not been appreciably depleted. This small region 
lies between the two convective regions: the one associated 
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to the inner He-burning shell and that associated to the 
H-rich envelope. In panel (b) it can be seen that the mix- 
ing due to the pulse-driven inner convective zone, allows 
a very small enrichment in 12 C and 16 of the regions just 
below the BCE. Panel (c) shows the maximum advance in- 
wards of the BCE, that is ultimately halted by the steep 
composition profiles that correspond to the products of 
the HBS. Finally, in panel (d), one can clearly see the ad- 
vance of the carbon profile, as a result of HeBS during 
the pulse, as well as the early recession of the BCE. In 
summary, even though overshooting has been taken into 
account, the amount of material dredged-up from the core 
to the convective envelope is very small. 

4.2. Models with oxygen-neon cores 

We will refer to models developing oxygen-neon cores as 
TP-SAGB stars. The upper and lower right panels of fig- 
ure[6]show the first thermal pulses of the 8 Mq model com- 
puted without overshooting and the corresponding ones 
for the 7 Mq model star computed with overshooting. 
These models have undergone extensive carbon burning 
and, consequently, have developed oxygen-neon cores. At 
the end of the carbon burning phase the core contracts 
while the helium burning shell is still active. As the fuel 
is exhausted in the innermost layers, this shell advances 
outward and approaches the base of the hydrogen-rich 
convective envelope. When the distance between the lat- 
ter and the HeBS is of about 10 -4 Mq the base of the 
hydrogen-rich convective envelope slightly recedes and the 
small radiative H-rich region is heated. Eventually the 
temperature in this region reaches the critical tempera- 
ture for the onset of the CNO-cycle and the HBS is re- 
activated. Therefore, the TP-SAGB phase develops in a 
very similar way to that of the models with carbon-oxygen 
cores. The details of the pulses of the 8 Mq model star 
computed without overshooting and of the 7Mq model 
in which overshooting was taken into account are shown 
in the upper and lower right panels of figure [71 Again, our 
calculations show that the enrichment in metals of the 
stellar surface is very small. 

Table [2] lists the surface abundances (by mass) of car- 
bon, nitrogen and oxygen, their relative abundance ratio 
(C:N:0) and the total metallicity (Z) of the model se- 
quences computed without overshooting (top section of 
table [J) and the corresponding ones when overshooting 
was taken into account (bottom section). As it can be 
seen, the total metallicity of the surface increases with in- 
creasing Mzams- The reason for this behavior is that more 
massive models have a lower degree of degeneracy in their 
cores and, therefore, the base of the convective envelope is 
able to penetrate deeper into the star and, consequently, 
is able to dredge-up larger amounts of material which has 
been processed by the different nuclear burning shells. 

Another important feature of the surface abundances 
presented in table [2] is that all models are carbon-rich, 
regardless of the adopted criterion for convection. For 
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Fig. 10. Masses of the carbon-oxygen cores after the cen- 
tral He burning phase for the models computed without 
overshooting (crosses and thick line) and with overshoot- 
ing (open circles and thin line). The mass coordinate for 
which carbon ignition occurs is also displayed for the mod- 
els computed without overshooting (triangles and thick 
line). Finally the mass thresholds at the ZAMS for car- 
bon burning to develop extensively are shown as vertical 
dotted lines. 

the models computed without overshooting (top section 
of table [2]) the relative abundance of nitrogen with re- 
spect to carbon first decreases, reaches a minimum at 
Mzams — 7 Mq and then increases again. The same hap- 
pens for the models computed with overshooting, but in 
this case the minimum is located at Mzams — 6 Mq . A 
similar behavior is obtained for the relative oxygen abun- 
dance with respect to carbon in the models computed 
without overshooting. For the case in which overshoot- 
ing has been taken into account the relative oxygen abun- 
dance shows the reverse behavior: it first increases up to 
a maximum at Mzams — 7 Mq and then decreases. The 
relative abundances of nitrogen and oxygen with respect 
to carbon, can be better understood if the models that de- 
velop carbon-oxygen or oxygen-neon cores are considered 
separately. For instance, the abundances of nitrogen and 
oxygen relative to carbon decrease as the Mzams increases 
for models with a carbon-oxygen core, and increase for the 
models which develop an oxygen-neon core. 

5. The properties of the cores 

Fig. [10] shows the mass of the carbon-oxygen core after 
the central He burning phase for the models computed 
with overshooting (open circles and thin line) and with- 
out overshooting (crosses and thick line) as a function of 
the mass at the ZAMS. The mass coordinate where carbon 
ignition initially occurs is also shown (triangles and thick 
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Table 2. Surface abundances (by mass) of carbon, oxygen and nitrogen at the end of our calculations for the model 
sequences computed without overshooting (top section) and with overshooting (bottom section). 
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line) for the models computed disregarding overshooting. 
Carbon ignition for the models computed with overshoot- 
ing is always central for this mass range. The ZAMS mass 
corresponding to the intersection of the thick lines sets a 
lower mass limit for single Z = stars to develop ONe de- 
generate cores. However, our calculations show that even 
though this is an actual lower limit to the mass thresh- 
old for the onset of carbon ignition, the real mass thresh- 
old for the formation of an ONe degenerate core is about 
2Mq larger, as stars of masses at the ZAMS between 
about 5.5 Mq and 7.8 Mq do not develop carbon burn- 
ing extensively enough in the carbon-oxygen core to let 
the core compositions be significantly changed. The dot- 
ted vertical lines correspond to these critical masses below 
which carbon is burnt only partially. For the case of the 
models computed with overshooting this mass threshold 
is located at Mzams — 6.0 Mq, whereas for the models in 
which overshooting was not taken into account the mass 
threshold is Mzams - 7.8 Mq. 



We have followed the evolution of all our model stars 
until the end of the carbon burning phase (if this is the 
case) including the early stages of the thermally pulsing 
AGB phase. Consequently, we have been able to deter- 
mine the chemical composition of the outer convective 
envelope and of the degenerate cores resulting from the 
evolution of massive intermediate mass stars, both for the 
case in which overshooting was considered and for the case 
in which it was disregarded. As it will be justified in sec- 
tion §6 with further detail — where we will discuss the 
possible evolutionary outcomes for this stellar mass range 
— the final destiny of the cores of these stars is not neces- 
sarily to form single white dwarfs. This only may occur if a 
mechanism exists able to produce extensive mass-loss and, 
hence, able to remove the convective envelope and leave 
behind a compact remnant. Such a mechanism could be, 
for instance, Roche lobe overflow due to the presence of 
a close companion, or, possibly, stellar winds powered by 
the (otherwise small) presence of metals in the envelope 
or, finally, by chromospherically driven winds (van Loon 
2005). Whatever this mechanism could be — if it indeed 
occurs — the mass and the chemical composition of the 
resulting compact remnants is also a by-product of our 
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Fig. 11. Composition profile of the degenerate CO core 
for the Mzams = 7Mq case without overshooting. 



calculations. We discuss these properties in the following 
subsections. 



5.1. Models computed without overshooting 

For the case in which no overshooting is adopted, primor- 
dial stars with masses at the ZAMS smaller than 7.8 Mq 
do not undergo extensive carbon burning. The pollution 
in metals of the stellar envelope is not large (see table [5]) . 
Hence, radiatively driven stellar winds during this stage 
are not likely to be strong enough to remove the outer en- 
velopes of these stars and produce carbon-oxygen white 
dwarfs, as it would be the case for solar metallicity stars 
of equivalent mass. Instead, the cores of these stars, due 
to the successive thermal pulses, could eventually grow 
up to the Chandrasekhar mass (M^h)- Thus, these stars 
are likely to undergo a supernova explosion. This possi- 
bility will be further explored and quantified in §6. If this 
is indeed the case, the chemical composition of the rem- 
nant core is of interest. As an example of the resulting 
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Fig. 12. Composition profile of the degenerate ONe core Fig. 13. Composition profile of the degenerate CO core 



for the Mzams = 8 Mq model without overshooting. 



for the Mzams = 5 Mq case with overshooting. 



degenerate core, figure [TT] shows the composition profile 
of the degenerate core resulting from a 7.0 Mq, computed 
without overshooting. 

According to our calculations, Z = stars in the range 
7.8 Mq <, Mzams < 9.0 Mq experience carbon burning 
in conditions of partial degeneracy if no overshooting is 
adopted. Thus, these stars form oxygen-neon cores sur- 
rounded by a carbon-oxygen degenerate buffer. For this 
mass range we also find that the enrichment in metals 
of the convective envelope is rather small and, hence, are 
prone to undergo accretion induced collapse (AIC) lead- 
ing also to a supernova explosion. The chemical compo- 
sition of such cores is also relevant for the particulars 
of their final evolutionary phases — see, for instance, 
Gutierrez et al. (1996) and Gutierrez et al. (2005). Fig. 
[T2l shows the internal chemical profile for the case of a 
Mzams = 8Mq primordial star. Finally, for the case 
in which no overshooting is adopted, primordial stars of 
masses Mzams i<2 9.5 Mq undergo extensive carbon burn- 
ing if no overshooting is taken into account, but the result- 
ing oxygen-neon cores are only weakly degenerate. Such 
stars would most likely become SNell directly. 

The top section of table [3] shows a summary of all this 
information. In particular, the mass of the CO core (sec- 
ond column) or the ONe core (third column) — depend- 
ing on whether the model has avoided carbon ignition or 
not — at the beginning of the thermally pulsing phase 
are displayed. Also shown are the masses of the remain- 
ing CO buffers on top of the ONe cores (AMco) for the 
models which have been able to ignite carbon (fourth col- 
umn), the carbon-to-oxygen ratio in the CO cores of the 
models which fail to ignite carbon (fifth column) and the 
neon-to-oxygen ratio for the models which ignite carbon 
(last column). As can be seen, for these evolutionary se- 
quences the mass of the carbon-oxygen buffer decreases 



as the mass of the ONe core increases and so does the 
neon-to-oxygen ratio. 

5.2. Models computed with overshooting 

The minimum mass at the ZAMS of a primordial star that 
allows for extensive carbon ignition and, hence, for the 
formation of an oxygen-neon core is, as previously men- 
tioned, 6.0 Mq. The convective envelopes of stars beyond 
this mass threshold are more metal rich than those com- 
puted without overshooting. Therefore, we expect that 
these stars could potentially be able to support radiatively 
driven winds and loose their envelopes, leading to the for- 
mation of carbon-oxygen white dwarfs. In this case the 
core composition of these stars is also of interest because 
the cooling speed depends on the carbon to oxygen ratio 
(Segretain et al. 1994). Figure Q2] shows an example of 
the resulting carbon-oxygen core for the Mzams = 5 Mq 
model sequence. 

Primordial stars with masses at the ZAMS between 
6.5 Mq and 8.5 Mq do ignite carbon at the center and 
may eventually become oxygen-neon white dwarfs, pro- 
vided that the convective envelope is removed by some 
mechanism during the thermally pulsing phase. We note 
however that in this case the metal content of the convec- 
tive envelope is sizeable — see table 2 — and, thus, this 
stars are prone to loose their convective envelopes as well. 
Fig. [14] shows the typical final chemical profiles of such 
an oxygen-neon white dwarf, which do not differ substan- 
tially from those of stars of solar metallicity. Finally, for 
stars with masses at the ZAMS larger than 8.5 Mq the hy- 
drogen burning shell does not resurrect and, thus, these 
stars do not have a thermally pulsing phase. Moreover, the 
mass of their oxygen-neon cores is larger than Mch and, 
consequently, their evolution continues all the way to more 
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Table 3. Characteristics of the remnant cores at the end of our calculations for the model sequences computed without 
overshooting (top section) and with overshooting (bottom section). 
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Fig. 14. Composition profile of the degenerate ONe core 
for the Mzams = 7 Mq model with overshooting. 

advanced stages of stellar evolution and they will most 
probably produce a supernova explosion. The bottom sec- 
tion of table [3] shows a summary of all this information. 
Finally, in figure [15] we present all this information graph- 
ically for both the case in which overshooting is adopted 
(bottom panel) or neglected (upper panel). In both panels 
the solid lines represent the total mass of the degenerate 
cores and the dotted lines the mass of the ONe cores. We 
have also separated and labeled the regions for which at 
the end of our calculations we obtain a CO core, an ONe 
core or directly a supernova. Note, however, that the final 
outcome of these stars still remains to be studied. This 
will be the goal of our next section. 

6. Mass loss and the final fate of primordial 
massive intermediate-mass stars 

In this section we will consider how the uncertainties that 
surround the evolution of primordial TP-(S)AGB stars 
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Fig. 15. Mass of the compact remnant cores by the end of 
our calculations for the cases with no overshooting (upper 
panel) and with overshooting (lower panel). The solid line 
shows the total mass of the degenerate core, while the 
dotted line is the mass of the ONe core — when applicable 
— without the CO buffer. 



can affect their expected final fate. We will show that, 
taking into account reasonable prescriptions for the mass 
loss rates during the TP-(S)AGB phase, the possibility 
that some of these stars end their lives as SNI1/2 — in- 
stead of white dwarfs, as it is usually assumed — remains 
open. These supernovae are produced by the explosion of 
a CO core in AGB stars (Arnett 1969, Nomoto 1976). We 
will also show that primordial stars that form ONe cores 
near the Chandrasekhar mass will possibly undergo an 
electron-capture induced collapse to a neutron star, very 
likely producing a type II-P supernova. Our study can be 
considered as a first approximation to the problem, since 
following the full evolution of these kind of stars until the 
very late phases of their evolution is clearly out of our cur- 
rent possibilities. Synthetic (S)AGB models like those of 
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Izzard et al. (2004) can help in this regard. Such an study 
is currently under way, and we postpone the discussion of 
such synthetic models to a forthcoming publication. 

The particular details of the evolution of stars of in- 
termediate mass from the end of the AGB phase, through 
a planetary nebula phase, until the formation of a white 
dwarf are not entirely well understood for solar metallicity 
stars. However, there is a general agreement on the overall 
picture. Strong winds are expected to remove most of the 
stellar envelopes before any significant growth of the cores 
due to shell burning processes is achieved. Therefore, the 
masses of the remnant cores after the main central burn- 
ing phases provide a good approximation to the final mass 
of the white dwarfs resulting from the evolution of these 
stars. 

However, the uncertainties associated to the mass-loss 
rates due to stellar winds affect the details of the late 
evolutionary stages of solar metallicity stars. For the case 
of primordial stars, these uncertainties are more dramatic 
because, as it will be shown below, the adopted mass-loss 
rates can change the final outcome of the stellar evolu- 
tionary sequences and let them end their lives either as a 
white dwarf or as a supernova. In principle, the mass-loss 
rates are expected to be smaller for stars of low metallic- 
ity. For instance, a common prescription to take this fact 
into account is to adopt the following expression: 



M(Z) = M(Z Q ) 



z. 



© 



0.5 



(2) 



where M(Zq) is the mass-loss rate for solar metallicity. 
Note, however, that this prescription is still rather un- 
certain and that it is the subject of a vivid debate (van 
Loon 2005). However, following closely the recent work of 
Ziljstra (2004), which is suggestive of a break-down of the 
mass-loss rates for dust-driven winds when Z < 0.1 Zq, 
we adopt it as an educated guess. Other recent works that 
have used this expression to scale mass loss with metallic- 
ity are Kudritzki (2000), Vanbeveren (2001) and Woosley 
et al. (2002). 

There are several prescriptions in the literature that 
provide M(Zq). One may use, for instance, the expression 
of Schroder & Cuntz (2005), whose authors claim it to 
yield reasonable results, even for the tip of the AGB: 



M sc (%) = Mr(Z ) 



cff 



4000 
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(3) 
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where T c ff is the effective temperature of the star, g is its 
surface gravity, <?q is the surface gravity of the Sun and 
A/r, is the mass loss rate of Reimers (1975): 

Mk(Zq) - -4 x 10~ 13 m ^ M yr- 1 (4) 



where t^r is a parameter such that 1/3 < r/p, < 3. 
Another common prescription is the expression proposed 
by Blocker (1995): 

Mb(Zq) = -4.83 x l O - 9 -^ T M R (Z ) M yr^ 1 (5) 

where Mpp is the actual mass of the considered star and, 
therefore, decreases as mass is lost. For solar metallicity 
stars this prescription yields consistently larger values dur- 
ing the TP-(S)AGB phase (Willson 2000, Gallart et al. 
2005) as compared to the rest of widely-used prescriptions 
found in the literature. We adopt it because it provides an 
upper bound for the expected mass-loss rate. 

When helium burning sets in a shell and during the 
thermally pulsing phase there is some dredge-up of metal- 
rich material and, consequently, the originally almost pri- 
mordial convective envelope becomes polluted in metals. 
Naturally, the metallicity of the convective envelope de- 
pends on the adopted prescription for convection, and, of 
course, on the evolutionary model, but it is never larger 
than Z ~ 2.8 x 10~ 3 . Taking into account Eq. (2), with the 
different prescriptions for M{Zq) given by Eq. (3), with 
ijn = 1.0, or by Eq. (5), a rough estimate of the minimum 
time required to loose the envelope through radiatively 
driven stellar winds can be obtained: 



tc 



AM, 



M 



(6) 



where AM onv is the mass of the convective envelope. In 
the bottom panel of figure ITrJI we have represented the log- 
arithm of the metallicity of the envelope, logZ onv , at the 
end of our calculations for the 5, 6, 7, 8 and 9 Mq stellar 
evolutionary sequences computed without overshooting 
(solid line), and for the 5,6,7 and 8 Mq model sequences 
computed with overshooting (dotted line). For the case 
of the 9 Mq model star computed with overshooting we 
did not obtain dredge-up and, thus, the envelope was not 
significantly enriched in metals, leading to Z cnv = 10 -11 . 
Additionally, the degenerate core of this model star has 
a mass very close to the Chandrasekhar mass and, thus, 
the most likely outcome is a supernova. The top panel 
of figure [16] represents the corresponding mass-loss rates, 
as obtained by adopting the prescription of Schroder & 
Cuntz (2005). We remark that these mass-loss rates could 
be taken as a lower limit, as they have been computed tak- 
ing into account the metallicity of the envelope at the end 
of our calculations. If the position of the BCE is, approx- 
imately, sustained and the envelope mass decreases, the 
freshly synthesized metals dredged-up from the core will 
appear less diluted in the stellar surface and hence will 
allow a certain increase in the mass-loss rates. On the 
other hand, the mass loss of primordial stars is actually 
hampered by the compactness and relatively high effective 
temperatures of these stars that result as a consequence of 
their low envelope metallicities. Since these stars are more 
compact the material lost by these stars must escape from 
a deeper gravitational well. Additionally, the high photo- 
spheric temperatures do not allow for the formation of 
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Fig. 16. Bottom panel: metallicity of the envelope of the 
selected models at the end of our calculations. Top panel: 
expected mass-loss rates given the envelope metallicities 
of the bottom panel when the Schroder & Cuntz (2005) 
mass-loss rates are adopted. Models computed without 
overshooting are represented using solid lines whereas 
models in which overshooting was adopted are shown with 
dotted lines. 

dust, even if refractory elements are present. Both facts 
will result into smaller mass-loss rates. 

Our evolutionary models yield the rate of growth of 
the degenerate core, M corc , as a consequence of the suc- 
cessive thermal pulses and also provide us with the mass 
of the degenerate core. Hence, the time needed to reach 
the Chandrasekhar mass would be: 



tc. 



M Ch - M c , 



(7) 



For instance, for the 7 Mq model star computed with 
overshooting the growth rate of the degenerate core during 
the thermally pulsing phase is ~ 2.0 x 10~ 7 Mq /yr. Since 
this model star has already an oxygen-neon core of mass 
MoNe — 1.2 Mq — see table[3] — it would only need about 
1.0 x 10 6 yr to reach the Chandrasekhar mass. Given the 
mass-loss rate shown in the top panel of ngureQl)] — about 
1.7 x 10~ 8 Mq yr -1 — the time needed to get rid of the 
envelope would be about 3.5 x 10 8 yr if the prescription 
of Schroder & Cuntz (2005) is adopted for the mass-loss 
rate. That is, t env is about 350 times larger than the time 
needed by the star to reach Mch- Therefore, one could 
reasonably expect that this star would end its life as a 
supernova. 

In this respect, it must be said that our work is just a 
first approach, as we have assumed that the rate of core 
growth is constant, that the third dredge-up is going to 
remain inefficient during the whole TP-(S)AGB phase, 
and that we have not included the effects of hot-bottom 



burning. A more realistic treatment of the problem would 
require to build up a synthetic code, as it has been done 
by Izzard et al. (2004), or Groenewegen & de Jong (1994) 
and references therein. However, there are some facts that 
support the interest of our results. First of all, the fact 
that the times required for the growth of the core up to 
Mch are several orders of magnitude longer than the times 
required to get rid of the envelopes. Second, the fact that 
even when disregarding the effect of metallicity, we still 
obtain the same qualitative behavior, although the values 
of the times required to eject the envelope decrease con- 
siderably. Third, that using other mass loss prescriptions, 
such as Reimers (1975) or Vassiliadis & Wood (1993), we 
also reproduce the same behavior. 

The procedure explained above can be repeated for 
all the evolutionary sequences presented here adopting ei- 
ther the prescriptions of Schroder & Cuntz (2005) and of 
Blocker (1995) for the mass-loss rates. The results can be 
found in table [4] The first column of this table displays 
the initial masses of the models analyzed in this work, 
the second column shows the time needed for the core 
of these stars to reach the Chandrasekhar mass, ^core 5 as 
obtained from Eq. (7), and using the core growth rates 
obtained from our evolutionary calculations. In particu- 
lar, the core growth rates have been obtained by dividing 
the mass growth of the core between the fifth and the last 
pulse computed in each model by the corresponding time 
interval. We have avoided the first pulses (or mini-pulses) 
in order to obtain average core growth rates during the 
whole thermally pulsing phase. With this procedure we 
obtain core growth rates ranging from 2 x 10~ 7 Mq yr -1 
to 5 x 10~ 7 Mq yr -1 . Because we have followed a limited 
number of pulses, these values could not correspond ex- 
actly to the average core growth rates during the whole 
thermally pulsing phase but provide us with a reasonable 
good approximation for the purposes of the present work. 
The time needed to get rid of the envelope according to 
the prescription of Schroder & Cuntz (2005), t^ v , is shown 
in the third column of tabled! whereas the time needed to 
loose the envelope according to the prescription of Blocker 
(1995), tg nv , is shown in the fourth column. As it has been 
done so far, the top section of table 0] shows the results 
for the models computed without considering overshoot- 
ing, whereas the results obtained taking overshooting into 
account are shown in the bottom section of this table. 

If we compare the times required for these stars to 
loose their envelopes according to the prescription of 
Schroder & Cuntz (2005) with those required by the cores 
to reach Mch, it can be seen that for all the models in 
table HlieS' ^ s between one and three orders of magnitude 
larger than t corc . Therefore, all these stars would likely 
end up their lives as supernovae, unless some additional 
mechanism enhancing the mass-loss rate played a signif- 
icant role in their evolution. Rotational induced mixing 
could be this extra mixing mechanism required for the ex- 
istence of strong stellar winds. Being completely devoid 
of metals, primordial stars would have no means of slow- 
ing down their initial rotation (Chiappini et al., 2006). 
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Table 4. Timescales (in years) associated to core growth 
and loss of the envelope for the case in which no over- 
shooting was taken into account (top section) and for the 
case in which overshooting was adopted (bottom section). 
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Furthermore, because they are more compact than higher 
metallicity stars of similar masses, higher rotational veloc- 
ities could be expected simply as a consequence of angu- 
lar momentum conservation during the initial collapsing 
stages. Hence primordial stars might experience a substan- 
tial amount of rotational mixing (Meynet et al., 2006). It 
is important to point that overshooting cannot account 
for this extra mixing, as rotational mixing proceeds also 
in radiative zones. In particular, it is important to note 
that Meynet et al. (2006) have followed the evolution of 
a 7Mq star, with Z = 10~ 5 , and have found a 1000-fold 
increase in its surface metallicity. Nevertheless, it is also 
worth noticing that the initial rotation velocity consid- 
ered in this work was near the critical rotation velocity. In 
any case, rotation is always difficult to introduce in a self- 
consistent way. Some groups assume that zero metallicity 
stars are fast rotators, but other advocate, on the basis 
of theoretical grounds, for much smaller rotation rates for 
primordial stars (Silk & Langer 2006). Moreover, Yoon & 
Langer (2005) point to another open question: the feed- 
back between rotation and stellar winds. These authors 
conclude that the rotation rates required to allow for the 
mixing of the collapsar can only be kept for weak stellar 
winds. Consequently, this point continues to be an open 
issue. 

The picture changes considerably when we consider the 
times required to loose the envelope when the prescription 
of Blocker (1995) is adopted. Not surprisingly, the mass- 
loss rates increase and, therefore, the times required to 
loose the stellar envelopes decrease. The results are the fol- 
lowing: the most massive models computed without over- 
shooting (the 8 and 9 Mq model stars) , and all the models 
computed with overshooting except the 8 Mq model are 
able to loose their envelopes before their cores reach Mch 
and, therefore, would probably end their lives as white 
dwarfs. The 8 Mq model computed with overshooting be- 
comes almost immediately a supernova, because the mass 
of its core (1.34 Mq) is already very close to Mch- The 
key factor that determines the long timescales required 
to loose the envelope for the models with masses between 
5 and 7Mq computed without overshooting is the low 



metallicities of the envelopes of these models. A compar- 
ison of the envelope timescales with the corresponding 
values of i coro allows us to conclude that these objects 
are likely to end their lives as SNeIl/2, as it has been 
previously found for all the cases computed adopting the 
Schroder & Cuntz (2005) prescription for the mass-loss 
rates. 

7. Summary and Discussion 

We have followed the evolution of massive intermediate- 
mass primordial stars from the zero age main sequence 
until the end of the carbon burning phase and the ther- 
mally pulsing (Super)-AGB phase (if this is the case). We 
have also studied the role of overshooting by computing 
two series of evolutionary sequences. In the first of these 
series overshooting was completely neglected. For this se- 
ries of evolutionary models we followed the evolution of 
model stars with masses ranging from 5Mq to 9Mq. In 
a second series of calculations overshooting was taken into 
account according to the prescription of Eldridge & Tout 
(2004) and we followed the evolution of primordial stars 
with masses ranging from 5 Mq to 8 Mq . This has al- 
lowed us to determine the masses of the degenerate cores 
resulting after the main central burning stages of primor- 
dial stars of these masses, and to study the mass limits 
for primordial stars to form CO white dwarfs, ONe white 
dwarfs or to undergo a supernova explosion. 

Our main results can be summarized as follows: if we 
consider the case in which overshooting was disregarded 
we have obtained a maximum mass for the formation 
of CO degenerate cores of ~ 7.8 Mq, and a maximum 
mass for the formation of an ONe cores of ~ 9.5 Mq. 
Stars more massive than 9.5 Mq form electron-degenerate 
cores with masses larger than the Chandrasckhar limit 
and, consequently, would directly undergo supernova ex- 
plosion through the core collapse mechanism. If a moder- 
ate amount of overshooting is adopted the mass threshold 
for the formation of a carbon-oxygen core, a neon-oxygen 
core or to undergo a supernova explosion is shifted to 
smaller masses by ~ 2Mq. Therefore, we get a maxi- 
mum ZAMS mass of 6.0 Mq for the formation of a CO 
degenerate core and a maximum value of 8.5 Mq for the 
formation of an ONe degenerate core. One would obtain 
a supernova for ZAMS masses larger than this value. For 
both sequences of models (those in which overshooting was 
taken into account and those in which it was disregarded) 
the resulting ONe cores appear to have a chemical strati- 
fication similar to that previously found by Garci'a-Berro 
et al. (1997) and Gil-Pons et al. (2005). Namely, an ONe 
electron-degenerate core surrounded by a relatively thick 
CO buffer whose absolute mass decreases as the mass of 
the degenerate ONe core increases. Typical masses for this 
CO buffer range from 0.04 to less than 0.01 Mq, depend- 
ing on the mass of the ONe core and on the choice of the 
prescription adopted for overshooting. The degeneracy in 
this CO buffer is also high. On top of the CO buffer there 
is a transition region in which helium is most abundant 
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and the degeneracy rapidly decreases and, on top of it, a 
H-rich convective envelope is found. 

The final outcome of these evolutionary sequences has 
been also analyzed in detail. Of course, this depends on 
both the rate of growth of the degenerate cores and on 
the adopted prescriptions for mass-loss. Moreover, since 
the mass-loss rates are most probably a function of the 
mctallicity of the convective envelope, the final outcome 
of these stars depends sensitively on this value and, thus, 
on the possible previous drcdge-up episodes. In partic- 
ular, for non-primordial stars such question would not 
represent a problem since stellar winds during the TP- 
(S)AGB stage are naturally pointed out as the mechanism 
by which these stars are deprived from their hydrogen-rich 
envelopes and, therefore, the formation of white dwarfs is 
immediately inferred. Such stellar winds can be easily jus- 
tified by the obvious presence of metals in the envelopes of 
non-primordial stars. However, only the lightest elements 
are originally present in Z = stars. Consequently, for the 
case of the massive intermediate mass stars that we have 
analyzed in this paper the situation is not so clear. Their 
evolution through core hydrogen , helium- and, eventu- 
ally, carbon-burning, and the dredge-up processes that 
accompany the TP-(S)AGB phase yield a modest enrich- 
ment in metals of the surface. Whether or not this en- 
richment is enough to enhance the mass-loss rates is un- 
clear. For instance, we have found that in the case of the 
7.0 Mq star computed disregarding overshooting the sur- 
face abundance by mass of metals is of the order of 10 -5 
after the 15 th pulse of the TP-AGB phase whereas in the 
case of the 9.0 Mq star the surface abundance of metals 
is already of the order of 10~ 4 by the end of the carbon- 
burning phase (Gil-Pons et al. 2005). The evolution dur- 
ing the TP-(S)AGB phase of the models computed with 
overshooting leads to more metal-enriched envelopes be- 
cause there exists a deeper penetration of the convective 
envelope to the regions of the star where core elements 
have been synthesized and, therefore, mass loss is more fa- 
vored. In spite of these metal enrichments we have shown 
that if the Schroder & Cuntz (2005) prescription (cor- 
rected for the effect of the metallicity of the envelope) is 
adopted for the mass-loss rates during the TP-(S)AGB 
phase, all the models studied here seem more likely to end 
up their lives as supernovae. The reason is that their cores 
are able to grow up to Mch in time scales between two and 
three orders of magnitudes shorter than those required for 
these stars to get rid of their envelopes, in agreement with 
the arguments of Zijlstra (2004). 

Taking the prescription of Blocker (1995) into account 
still allows the formation of SNI1/2 as a consequence of 
the evolution of ZAMS objects between 5 and 7Mq com- 
puted without overshooting. This prescription is very of- 
ten considered to produce very high mass-loss rates and, 
therefore, the fact that even under these conditions the 
possibility for the formation of SNI1/2 remains open is 
an encouraging result which deserves further scrutiny. In 
any case, it must be kept in mind that the lack of a first 
principle theory of mass-loss results in a considerable com- 



plication for these studies, even in the case of solar metal- 
licity stars. The conclusions on the subject must rely on 
detailed observations of open clusters or detached binary 
systems in which one of the members is a white dwarf — 
see Weidemann (2000) for a thorough review — and any 
extrapolation to primordial objects might be still some- 
how speculative. 
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